Order from disorder: Quantum spin gap in magnon spectra of LaTiOs 
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A theory of the anisotropic superexchange and low energy spin excitations in a Mott insulator with 
t2 S orbital degeneracy is presented. We observe that the spin-orbit coupling induces frustrating Ising- 
like anisotropy terms in the spin Hamiltonian, which invalidate noninteracting spin wave theory. 
The frustration of classical states is resolved by an order from disorder mechanism, which selects 
a particular direction of the staggered moment and generates a quantum spin gap. The theory 
explains well the observed magnon gaps in LaTi03. As a test case, a specific prediction is made on 
the splitting of magnon branches at certain momentum directions. 
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The recent renaissance in the study of transition metal 
oxides has emphasized the important role being played by 
the orbital degeneracy inherent to perovskite lattices (see 
for review |]||). First of all, the type of spin structure 
and the character of spin excitations crucially depend 
on the orientation of occupied orbitals Hf§. Second, 
the excitations in the orbital sector get coupled to the 
other degrees of freedom (electronic, lattice, spin) and 
might therefore strongly modify their excitation spec- 
trum. Some examples are the anomalous magnon soften- 
ing || and incoherent charge transport |J due to low- 
energy orbital fluctuations in ferromagnetic manganites, 
and the formation of orbital polarons [Q . 

The calculation of the spin excitation spectrum in sys- 
tems with orbital degeneracy is somewhat involved even 
at the half filled, insulating limit. That is due to the 
peculiar frustrations of superexchange interactions 
which lead to the infrared divergences when a linear spin 
wave theory is applied. In systems with e g orbital degen- 
eracy the frustration is resolved by a specific, directional 
orbital ordering thus breaking the cubic symmetry of spin 
exchange bonds Q . The transition metal oxides with tn g 
orbitals exhibit different and more challenging phenom- 
ena. This occurs due to the relative weakness of the 
Jahn- Teller coupling in this case, and due to the higher 
degeneracy and additional symmetry of t2 g orbitals [fl0[ . 
As a result, the orbitals may form the novel, coherent 
orbital-liquid state stabilized by quantum effects, as ob- 
served in spin the S — 1/2 Mott insulator LaTiOs 

LaTiOs shows a spin order of G-type with magnon 
excitations characteristic for a simple cubic lattice. De- 
spite three dimensionality of magnon spectra, the spin 
reduction is unusually large (much larger than in two di- 
mensional cuprates), which has been explained in terms 
of fluctuating orbital state in this material Jll]jl(| • The 
present paper concerns with the origin of the small spin 
gap observed in LaTi03. This would not be an issue in 
a conventional case with static orbital order: The lat- 
ter lowers a symmetry of the crystal and induces a spin 
anisotropy via the spin-orbit coupling, resulting naturally 
in a classical gap in the magnon spectra. Orbital order 



is however not observed in LaTiOs, and the way how the 
staggered moment chooses one out of three equivalent 
cubic axes as the easy one is not that obvious. Indeed, it 
is shown below that the anisotropic Hamiltonian induced 
by spin-orbit coupling has a perfect cubic symmetry in 
the orbital liquid state, and a conventional spin wave 
theory gives in fact no magnon gap. We argue that the 
magnon gap and the selection of the easy magnetization 
axis in LaTiOs are nontrivial effects of quantum origin 
generated by the order from disorder mechanism. 

We begin with the derivation of the effective spin 
Hamiltonian in i 2 g systems. Quite generally, it consists 
of the isotropic superexchange of Heisenbcrg form, and 
the anisotropic spin exchange Hamiltonian. The latter, 
which is of our present interest, results from higher order 
processes involving both the spin-orbit coupling 
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and the isotropic superexchange. Here the operator 1; is 
the effective angular momentum (I = 1) of t2 g level, and 
S being a spin one-half of the Ti 3+ ion. The isotropic 
superexchange in orbitally degenerate systems strongly 
depends on the orbital structure. In general it can be 
written as: 
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where the orbital operators J^- and depend on 

bond directions j(= a,b,c). In a t2 g system like the 
titanates they are given by the following expressions: 
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The coefficients r\ — 1/(1 — 377) , r 2 = 1/(1 — 77) and 
rs = 1/(1 + 2rf) originate from the Hund's splitting of the 
excited t\ g multiplet via 77 = Jh/U. The operators 
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FIG. 1. On every bond of cubic crystal, two out of three 
ti g orbitals are equally involved in the superexchange and 
may resonate. The same two orbitals select also a particular 
component of angular momentum. 



and nf' can conveniently be represented in terms 
of spinless fermions a i; a corresponding to t 2g levels 
of yz, xz, xy symmetry, respectively. (This notation is 
motivated by the fact that each t 2g orbital is orthogonal 
to one of the cubic axes a,b,c). Namely, 



,(7) 
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for the pair along the c axis. Similar expressions are ob- 
tained for the exchange bonds along the axes a and b, 
by replacing orbital fermions (a, b) in Eqs.(5) by (6, c) 
and (c, a) pairs, respectively. Since the parameter 77 = 
Jh/U ~ 0.16 only we will neglect the Hund's cou- 
pling corrections. This results in a simple expressions as 
obtained in jig ]: 

if 7 c) = n ia n ja + n lb n jb + alb^a^ + l>]u,a"b,. (6) 
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The Eq.(7) is actually unessential for our purpose. 
Summed over bonds, this term gives the energy of the 
classical Neel state. 

To complete the definitions, we express the angular 
momentum operator in Eq.(l) as follows: 

l\ = i{a\b t - ftjoi), If = i(b\a - c\b,),l v t = i{a\a - cfoi), 

(8) 

and notice finally that the fermions must satisfy a local 
constraint, rii a + ria + rii c = 1. 

To proceed further we need to specify the orbital state. 
The two most important features of the interactions in 
the present model can be observed from Eqs.(2,6): z)The 
classical Neel state (consider (S,Sj) = —1/4) is infinitely 
degenerate in the orbital sector, thus strong spin fluc- 
tuations must be present in the ground state to pro- 
vide a dynamical splitting of orbital levels. n)Special 



to t 2g systems, every bond is represented by a particular 
pair of equivalent orbitals (see Fig.l and Eq.(6)), which 
may form dynamical orbital singlet or triplet bonds. (In 
fact, the orbital operators may be represented in a 
SU(2) symmetric form [[To||). This additional symmetry 
brings up an intrinsic quantum dynamics into t 2g orbital 
physics, contrasting it to a classical behaviour of e g sys- 
tems. As described in JTo[ ] , a composite spin-orbital ex- 
citation (an analogy to a SU(4) excitation jl3|]), which 
consists of a spin fluctuation accompanied by the for- 
mation of dynamical orbital bonds, best optimizes the 
superexchange energy. In that orbital liquid picture cu- 
bic symmetry remains unbroken. Angular momentum 1 
is however fully quenched in the ground state, and it's 
fluctuation spectrum extends up to characteristic ener- 
gies of about W orb ~ 4:t 2 /U. Last important remark 
concerns with spin-orbital separation: While the large 
energy (~ W or b) behaviour of the model Eq.(2) is gov- 
erned by a coupled spin-orbital dynamics, the underlying 
weak magnetic order results in a separation of low energy 
spin excitations forming magnons on scale J. Here J is a 
pure spin exchange coupling, J(SjSj), which is obtained 
by averaging Eq.(6) over fast orbital dynamics. It has 
been found that indeed J <§; W or b JTofl , and this justi- 
fies the adiabatic approximation used below to derive an 
effective spin anisotropy Hamiltonian. 

By symmetry, the anisotropic pairwise interactions for 
spins one-half in a cubic crystal must have the form 
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This interaction follows in the present model from spin- 
orbit corrections to the superexchange, shown in Fig. 2. 
Only two of them, diagrams a) and c) do actually con- 
tribute to the anisotropy constant A. Consider the pair 
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Here (• • -) or b implies the average over orbital fluctua- 
tions. The basic assumption in the above derivation 
is that these fluctuations are fast enough to enable us 
to integrate them out as far as we are concerned with 
the low energy, magnon like dynamics of the spin op- 
erators. Next, it is important to observe that both l z 
(see Eq.(8)) and jj. operate on the same orbital dou- 



blet (a, 6). One therefore obtains A{ c ^ > A^' — A y c> . For 
the pairs (ij) \\ a and (ij) \\ b, on the other hand, (6, c) 
and (c, a) doublets are more active, resulting in ASfSJ 
and ASf Sj interactions, respectively. The physical pic- 
ture is that the anisotropic spin interaction arises here 
as an indirect coupling via fluctuations of the angular 
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FIG. 2. Lowest order spin-orbit corrections to the spin ex- 
change Hamiltonian for a pair (ij). Wavy lines represent the 
isotropic superexchange in Eq.(2), solid lines are the orbiton 
on-site Green's functions. The spin-orbit coupling in Eq.(l) 
is denoted by crosses. 

momentum of t 2g level. In every direction one partic- 
ular component among l a in Eq.(8) is special, having 
more dispersive fluctuations in that direction, therefore 
the structure of Eq.(9) follows. 

Eq.(lO) can be expressed in terms of an on-site orbiton 
Green's functions which we evaluate using constant den- 
sity of states within the orbiton bands of width W or b, as 
obtained in the large N approximation p0| . This gives 

A = f{At 2 /U){X/W orb ) 2 , (11) 

with the numerical factor / ~ 3. We may notice, that 
W or b in this expression actually plays the role of a static 
splitting of levels in a conventional, orbitally ordered 
case. The estimation in Ref. |l(| yields for the Heisen- 
berg exchange constant J ~ 0.16(4i 2 /£7) and the orbiton 
bandwidth W or b — l-6(4£ 2 /£7), thus the energy scale 
At 2 /U ~ 100 meV follows using J eX r, = 15.5 ± 1 meV 
in LaTi0 3 @. With A = 19 meV |j, we then estimate 
A ~ 4.2 meV for cubic anisotropy constant in LaTi03. 

We point out here a crucial difference between the 
present case and a canonical theory of the anisotropic su- 
perexchange interactions in orbitally nondegenerate sys- 
tems like cuprates. In e g systems, one has to treat both 
symmetric and antisymmetric (Dzyaloshinskii-Moriya 
(DM) jlqQ parts of the anisotropic superexchange on 
equal footing because of the intimate relationship be- 
tween them (161 ) since both are due to the finite devia- 
tions (tilting of octahedra etc.) from the ideal structure. 
In the present t 2g case these interactions are, however, 
decoupled: while the DM interaction vanishes by cubic 
symmetry, a symmetric superexchange anisotropy is fully 
operating. The key point here is that the t 2g state man- 
ifold has an intrinsic angular momentum. 

Now we consider how the interaction Eq. (9) will affect 
magnon excitations. Since spin symmetry is lowered to 
the discrete (cubic) one, a magnon gap is expected. A 
linear spin wave theory fails however to obtain it, because 
Eq.(9) acquires a rotational symmetry in the limit of clas- 
sical spins. This results in an infinite degeneracy of clas- 
sical states, and an accidental pseudo Goldstone mode 
appears, which is however not a symmetry property of 
the original quantum model Eq.(9). Villain's order from 
disorder mechanism |l7| comes into play at this point, 
selecting a particular classical state such that the fluctu- 



ations about this state give the largest energy gain, and 
opening also a magnon gap. To explore this point explic- 
itly, we calculate quantum corrections to the ground state 
energy as a function of the angle 9 between c axis and 
the staggered moment. Assuming the latter is perpen- 
dicular to b axis, we rotate globally a spin quantization 
axes as S z -> S z cos 9-S x sin 9, S x -► S x cos 9 + S z sin 9, 
and observe that the magnon excitation spectrum has an 
explicit ^-dependence: 

lu p (9) = zJS{Y p (X p + 2an zx sin 2 9)} 1 ' 2 , (12) 
Y p = 1 - 7 P + o(l - c y ) , (13) 
A p = l + 7 P + a(l + c x ) . (14) 

Here the ratio a = A/ 3 J quantifies the deviation from 
the Heisenberg limit, and r\ zx = (c z — c x )/2 , 7 P = (c x + 
c y + c z )/3 , where c Q = cosp a , and z = 6. Calculating 
the zero point magnon energy, one obtains an effective 
potential for the staggered moment, which at small 9 
reads as a harmonic one: E(9) = const + K e f j S 2 9 2 , with 
an effective "spring" constant K e ff — \zJan = An. The 
parameter re is given by 

k=^5>-(V*p) 1/2 - (15) 
p 

At small a <C 1, in which case we are interested in, re ~ 
30^/85, and K eff ~ (R/8S){A 2 / J), with the numerical 
factor 

^ = E( 1 +^p)( c '-^p)/( 1 -^) 3/2 -0.44. (16) 
p 

The physical meaning of the above calculation is that 
zero point quantum fluctuations, generated by the in- 
teraction in Eq.(9), are enhanced when the staggered 
moment stays about a symmetric position (one of three 
cubic axes), and this leads to the formation of the en- 
ergy profile of cubic symmetry. A breaking of this dis- 
crete symmetry results then in the magnon gap, which 
should be about y/K e ff/M with 1/M ~ zJ in the 
harmonic approximation. More quantitatively, the po- 
tential E(9) can be associated with an effective uni- 
axial anisotropy term, K e f fY2(ij) c generated in 
the symmetry broken phase. Therefore, one finds the 
magnon gap A = z JS^j2K e f ~t /3 J, which is a linear func- 
tion of the anisotropic interaction A: 

A = zJSV2ok = AV3RS . (17) 

Now we derive this result using a different approach. 
Namely, we write the equations of motion, uS° = [S° , H] 
with H being the Heisenberg exchange plus H an i in 
Eq.(9), and account for the interaction between magnons 
and local fluctuations of the staggered moment 5S Z , gen- 
erated by Hani-, within the following approximation: 
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Here 7 = a,b,c specifies the bond direction (ij), and the 
bond variables are defined as x^ ~ > A*- 7 ** = 

— (b\bj) , where 6, is the Holstein-Primakoff boson oper- 
ator. This linearization procedure leads to the magnon 
excitation spectrum 



= z JS{(Y p + aK)(X p + an)} 1/2 



(19) 



where n = - n^ c \ and kW = ( x h) + A^)/S. The 
key observation here is that the bond variables along the 
c axis differ from those in directions perpendicular to the 
staggered moment. Indeed, calculating bond variables 
and A^ by using the spectrum in Eq.(19), one finds 
k to be positive, given by the sclfconsistcnt equation: 



« = ^E VzxiYp + aK^^/iXp + arc) 
p 



1/2 



(20) 



This is expected since the fluctuations originate mostly 
from a and b bonds in Eq.(9), once we have chosen the c 
axis as an easy one. Thus rc^ > rc^ follows. 

Consider again the weak anisotropy case, a -C 1. From 
equations (19) and (13,14) we read off the magnon gap 
A = zJS\/2an. Noticing also that Eq.(20) reproduces 
in this limit precisely the same parameter k defined in 
Eq.(15), we confirm finally our previous result, given by 
Eq.(17), for the magnon gap. Superior to the previous 
case, the latter solution gives not only the gap value, but 
also the full magnon spectrum Eq.(19), which we will 
discuss soon. 

One may wonder about the behaviour of Eqs.(20,19) in 
another extreme limit, o> 1, where one is left with only 
the cubic term. Actually, the model defined in Eq.(9) 
is of certain interest by it's own jul. For a — > 00, 
one finds that the solution of Eq.(20) saturates at the 
value Koo ~ (2/n 2 S)rs, where r$ = ln(47r\/S) — 1- 
The excitation spectrum in this limit is two-dimensional, 
ujp = 2AS{(1 + — Cy)(l + Hoc + Cx)} 1 / 2 , and it shows 
again the gap A = (4 / V) Ay V s S of order A. The spin 
reduction due to the frustrations is found to be non- 
singular for finite spin values. Although it is logarith- 
mically diverging in the limit of large classical spins, 
SS Z — ► Jj- ln(47iV5) , the ratio SS Z /S vanishes in this 
limit. Therefore, we believe that long range order in the 
"cubic" model with nearest-neighbor interactions is well 
defined for both one-half and large spins. It should be 
noticed, however, that the approximation used here is 
less accurate in the "cubic" limit, and we anticipate sub- 
stantial incoherent features in the excitation spectrum, 
which is expected to show a pseudogap of about A and 
a full bandwidth u> max ~ 4^45". 

Now we turn to the realistic case, LaTi03, and dis- 
cuss the numbers. With A ~ 4.2 meV estimated above, 



Eq.(17) gives A ~ 3.4 meV, which is in surprisingly good 
agreement with A exp = 3.3 ± 0.3 meV |11 . We believe 



therefore that the quantum gap induced by the cubic 
anisotropy is the dominant contribution to the magnon 
gap in LaTi03 [|l9|. This can actually be tested ex- 
perimentally. The point is that the spectrum uj p in 
Eq.(19) contains full information of the interaction in 
Eq.(9) on large energy scales, and it's peculiar structure 
may manifest itself in a neutron scattering experiment. 
Namely, two magnon branches of a commensurate anti- 
ferromagnet are related to Eq.(19) as f2i(p) = lo p and 
£7 2 (p) = Wp+Q, where Q being the Neel vector. Because 
of the presence of S x and S y components in Eq.(9), the 
two branches are in fact degenerate only in the planes 
with p x = ip y , otherwise a splitting of these two modes 
is expected at large momenta. The splitting is largest 
along [0, 7r,p z ] and equivalent lines and reaches a maxi- 
mal value 2A ~ 8 meV at the p z = it/2 point. This may 
enable one to measure directly the value of the cubic 
anisotropy constant A, and quantify further the orbital 
liquid state in a such puzzling material like LaTiC>3. 

We would like to emphasize that a small value of the 
spin gap in LaTi03 is very difficult to explain from the 
traditional point of view (see also the discussion in jll| ) . 
Indeed, large static orbital splitting of order of 150-200 
meV would be required in order to suppress spin-space 
anisotropy and hence to explain the small magnon gaps 
as observed. However, such large splittings are very dif- 
ficult to reconcile with the fact that neither orbital or- 
der nor any traces of related structural transitions have 
been detected. The orbital liquid concept naturally re- 
solves this apparent conflict: The large orbital energy 
scale is generated in this picture by electronic correla- 
tions, and angular momentum is quenched dynamically, 
without static symmetry breaking. 

To conclude, a theory of the anisotropic superexchange 
of tig electrons in titanites is presented. This interac- 
tion introduces nontrivial frustration in low energy spin 
states. We have identified a mechanism which selects the 
easy magnetic axis in cubic crystal with disordered or- 
bital states. The spin gap in LaTiOs is most likely of 
quantum origin, being thus as unique as the orbital state 
in this "simple" Mott insulator is. An experimental test 
of the theory is suggested. 
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